Comparison of the phase diagram of the half-filled layered organic superconductors 
with the phase diagram of the RVB theory of the Hubbard— Heisenberg model 
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We present an resonating valence bond (RVB) theory of superconductivity for the Hubbard- 
Heisenberg model on an anisotropic triangular lattice. We show that these calculations are consistent 
with the observed phase diagram of the half-filled layered organic superconductors, such as the 
/3, fi\ k and A phases of (BEDT-TTF) 2 X [bis(ethylenedithio)tetrathiafulvalene] and (BETS) 2 X 
[bis(ethylenedithio)tetraselenafulvalene]. We find a first order transition from a Mott insulator to 
a d x 2_ y 2 superconductor with a small superfluid stiffness and a pseudogap with d x 2_ y 2 symmetry. 
The Mott-Hubbard transition can be driven either by increasing the on-site Coulomb repulsion, U, 
or by changing the anisotropy of the two hopping integrals, t' /t. Our results suggest that the ratio 
t'/t plays an important role in determining the phase diagram of the organic superconductors. 
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Describing strongly correlated electronic systems is one 
of the outstanding challenges of theoretical physics. In 
particular one would like to understand if different model 
materials embody the same underlying physics. The sim- 
ilarities between the cuprates and the layered organics 
superconductors 0, suggest that similar physics may 
be realised in both classes of materials. A powerful ap- 
proach to chemically complex materials, such as organic 
superconductors, is to define minimal models J3|, which 
can then be treated at various levels of approximation 
0, HI- In this Letter we take such an approach. We 
argue that the observed phase diagram of the half-filled 
layered organic superconductors (^LOS) is well described 
by the RVB theory of the Hubbard-Heisenberg model. 
Our theory reproduces the first order Mott transition and 
predicts d^2_ y 2 superconductivity 0, a small superfluid 
stiffness and a pseudogap 0. 

Layered organic superconductors form several crystal 
structures, some of which, such as the /3, (3', K and A 
phases are strongly dimerised, others, e.g., the a, f3" and 
9 phases are not. The chemical composition of these ma- 
terials is D2X where D is an organic donor molecule, for 
example BEDT-TTF (ET) or BETS, and X is an anion. 
Crystals consist of alternating layers of donor molecules 
and anions In both the dimerised and undimerised 
salts the anion accepts one electron from a pair of donor 
molecules which leads, at the level of band structure, to 
an insulating anionic layer and a metallic donor layer. 
Quantum chemistry suggests that the band structure of 
the undimerised materials is well described by treating 
each donor molecule as a site in a (quarter-filled) tight- 
binding model 01 • I n the dimerised materials the intra- 
dimer hopping integral is large enough that the band 
structure can be described by a half-filled tight-binding 
model with each site representing a dimer 0,0). 

|LOS display insulating, metallic, superconducting, 



being proposed. Both phononic and non-phononic pair- 
ing mechanisms have previously been considered 000- 
The superfluid stiffness is much smaller than is pre- 
dicted by BCS theory but is too large for fluctuations in 
the phase of the order parameter to be important • 

Fig. n shows the phase diagram of k-(ET)2A as a 
function of pressure (both hydrostatic and 'chemical') 
and temperature. Other ^LOS have similar phase di- 
agrams |9(. A simple explanation of this phase diagram 
is as follows 0, 0: there is a strong on-site Coulomb 
repulsion, U, which causes the ambient pressure (Mott) 
insulating state. The application of hydrostatic pressure 
or varying the anion (often thought of as applying 'chem- 
ical pressure') reduces U /W, where W is the bandwidth, 
and leads to a superconducting state caused by strong 
electronic correlations. The bad metal phase is due to 
somewhat localised electrons as one crosses over from 
the Fermi liquid to the Mott insulator (which does not 
require a phase transition in these materials |l5lll6| ). 

It has been argued that the Hubbard model on an 
anisotropic triangular lattice is a minimal model for the 
layered organic superconductors • A Dynamical Mean 
Field Theory (DMFT) of the Hubbard model on a hyper- 
cubic lattice gives a good quantitative description of the 
competition between the Mott insulator, the bad metal 
and the Fermi liquid |l6l Il9| . However, a mean field 
treatment of the positive U Hubbard model will not cor- 
rectly describe the materials as it neglects important spin 
correlations which arise from superexchange. We there- 
fore consider the Hubbard-Heisenberg model, which can 
be derived [2$ fr om the Hubbard model in the limit of 
large, but finite, U. The Hamiltonian is 

irr {ij}<r (ij)cr 

+J S» • S 3 +J'J2s i -S j + Uj2 n ir hu (1) 



'bad metallic' and (possibly) pseudogap 8] phases. The 
nature of the superconducting state in ^LOS is contro- 
versial 0] : the pairing is thought to be singlet @ , but ex- 
periments have lead to both s-wave and d-wave scenarios 



where c.-t 



annihilates (creates) an electron on site i with 
spin tr, Si is the Heisenberg spin operator, hi a is the num- 
ber operator, and {ij} and (ij) indicate sums over nearest 
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FIG. 1: (Color online.) The pressure-temperature phase dia- 
gram of k-(ET) 2 X. Top: data from *H NMR and AC suscep- 
tibility (dark blue diamonds show the transition from a non- 
magnetic state to an antiferromagnetically ordered state, light 
blue diamonds show the metal-insulator transition 12|), mag- 
netisation (pink pluses E3l)> thermal expansion (filled purple 
circles and resistivity (red squares {filled indicates a 

first order Mott transition, empty indicates a crossover from 
insulating to metallic behaviourslJl5|. filled gree n triangles 
flrjl . grey stars 0] an d open purple diamond |l8j 1. We have 
offset the data to allow for the effect of 'chemical pressure'. 
P = corresponds to ambient pressure for X =Cu[N(CN)2]Cl 
|I^llall4lllMll6|| . The 'chemical' pressure is indicated by the 
arrows on the abscissa. H 8 -NCS X =Cu(NCS) 2 [3, 
and H 8 -Br => X =Cu[N(CN) 2 ]Br EHH. D 8 -Br indicates 
the effective chemical pressure of X =Cu[N(CN)2]Br with 
the ET molecule fully deuterated fl^l . Bottom: a schematic 
version of the same diagram is shown. The glassy transition 
between the ethylene liquid and frozen in disorder phases re- 
sults from conformational disorder in the organic molecule 
6, 14]. 



neighbours and next nearest neighbours across one diag- 
onal only 3 respectively. In principle J = At 2 /U and 
ji _ 4^/2 jjj |- G leading order due to superexchange. How- 
ever, our mean-field treatment will not correctly describe 



the renormalisation of the bare parameters. Therefore 
we treat t 7 t' ', J, J' and U as independent parameters. 
To reduce our parameter space we choose J = t/3 and 
,/' = t' 2 /3t jU, which correspond roughly to the values 
of J and J' extracted from experiments on the insulating 
phase of the layered organics 0, |2^] . Thus in the calcula- 
tions presented below we only vary two parameters: t' jt 
and U/t. 

Our treatment of the Hubbard-Heisenberg Hamilto- 
nian Q is based on Anderson's RVB theory Al- 
though, the RVB wavefunction is a poor approximation 
for the Heisenberg model on a square lattice, it was re- 
cently shown that it is a good trial wavefunction for 
some frustrated Heisenberg models '24\. These models 
are closely related to ours as many of the ^LOS are 
expected 0, to have J'/ J = (t' /t) 2 in the relevant 
range. Further evidence that the RVB theory is a much 
better theory for the triangular lattice than it is for the 
square lattice comes from the critical value of the on site 
Coulomb repulsion, U c , at which the Mott transition oc- 
curs. The square lattice is insulating for arbitrarily small 
values of U, whereas we find that the RVB theory gives 
U c ~ 10. 3i. On the isotropic triangular lattice exact di- 
agonalisation of finite lattices gives U c = 12t and we 
find that U c ~ 12.4t in the RVB theory. 

Anderson's RVB state, \RVB) 1 is given by perform- 
ing a Gutzwiller projection, Pq = Y]Al — an^n^), on 
the BCS wavefunction, \BCS), i.e., \RVB) = P G \BCS). 
Here a is a variational parameter which controls the frac- 
tion of doubly occupied sites, d. A detailed analysis of 
the RVB theory of the Hubbard-Heisenberg model on 
the square lattice was reported by Gan et al. |2fjj ]. 

Following the spirit of Ref. we make the Gutzwiller 
approximation [2f|, viz., (c^Cj^rvb = 0i(4o-Cj>) bcs 
and (S t -Sj) RVB = gs(Si-Sj)Bcs where (O),/, = (ip\<D\ip). 
Counting arguments show that [2fjL |2S| at half filling 
g t = 8(1 - 2d)d and g s = 4(1 - 2d) 2 . The Gutzwiller 
approximation has several advantages: its simplicity al- 
lows some analytic progress to be made and allows one 
to consider infinite systems. However, the Gutzwiller ap- 
proximation suppresses spin and charge fluctuations in 
the Hubbard model [isf. We have already sidestepped 
this problem somewhat by explicitly including the spin 
exchange terms in the Hubbard-Heisenberg model. Our 
theory produces a Mott insulating state that is a spin liq- 
uid rather than the antiferromagnetic insulating state ob- 
served in most (but not all j^) |LOS, however general- 
isation of \RVB) to allow for antiferromagnetism should 
not significantly alter the phase diagram. Clearly an im- 
portant test will be to project \RVB) onto the results of 
exact diagonalisation of finite systems for the Hubbard 
model on the anisotropic triangular lattice. 

Making the Hartree-Fock-Gorkov approximation leads 
to two coupled gap equations, Ak = — X)k'^k- k '^T 
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FIG. 2: (Color online.) g t ~ Z as a function of the on- 
site Coulomb repulsion U for various levels of values of t'/t, 
throughout J = t/3 and J' = t' 2 /St |2ll |. Here we plot 
i' = 0.7£ (dashed), t' = 0.3* (dot dashed) and i' = (solid). 
For each value of t'/t we find a first order Mott-Hubbard tran- 
sition at some critical value of the on-site Coulomb repulsion, 
U c , from a superconducting state (gt =fi 0) to an insulating 
state (gt = 0). The inset shows the phase diagram for the 
model. Note that the fact that U c varies with t'/t shows that 
increasing t'/t can drive the Mott-Hubbard transition. 

— ^gs[J(cosk x + cosky) + J' cos(k x + k y )], £ k = p, — 
g t 2[t(cosk x + cosky) + t' cos(k x + k y )] and = 
\/[^k + Xk] 2 + |A k | 2 . d is minimised variationally and 
the renormalised chemical potential, /j, is chosen to en- 
sure half-filling 26]. The two mean- fields are a Hartree- 
Fockterm, Xk = J2k' Vk-k' (^k't^k't) BCS , and an anoma- 
lous term, A k = J2w Vk-k' (ck'TC-k'i)scs, where c k(T is 
the Fourier transform of Ci a . 

We solve the coupled gap equations self consistently in 
reciprocal space on a 120 x 120 mesh. We do not enforce 
any symmetry constraints on the order parameters and 
we find that it has d x 2_ y 2 symmetry, this is the pairing 
symmetry most compatible with a range of experiments 
on the layered organics |f|. 

For the Hubbard model gt = Z, the quasiparticlc 
weight (the factor by which many-body effects reduce 
the bandwidth and Drude weight and enhance the ef- 
fective mass, m*). In Fig. 0we plot g t against U for 
several values of t'/t. For all values of t'/t we find that 
at some critical value, U c , there is a first order transition 
from a superconducting state (gt ^ 0) to an insulating 
state (gt = 0). This is consistent with the first order 
superconductor-insulator transition observed experimen- 
tally iniLOSlUGi. 

Note that near the Mott transition 
gt is reduced (and hence m* is enhanced by a factor of 3 
or 4), consistent with the large effective mass seen in the 
layered organic superconductors close to the insulating 
state Q. Previous weak coupling approaches 0] do not 
capture this large mass renormalisation. gt is closely re- 
lated to the reduction in the Drude weight due to strong 
correlations. Fig. 2 is quantitatively similar to to exact 



FIG. 3: (Color online.) The mean gap, |Ak|, and the mean 

superconducting order parameter, | Ak], (both in units of t) 
as functions of the on-site Coulomb repulsion U for various 
values of t'/t, throughout J = t/3 and J' = t' 2 /3t 2lj. A k = 
<?t A k and the bar indicates averaging over the Brillouin zone 

(i.e., ~]Z = Ek/k/Ek)- Here we p lot t' = 0.7* (blue: |A k | 

solid; | Ak | dashed), t' = 0.3t (red: |Ak| double dot dashed; 

|Ak| long dashed) and t' = (black: |Ak| dot double dashed; 

|Ak| dot dashed). The inset shows |Ak|, (dashed) and |Ak|, 
(solid) against t'/t with J = t/3, J' = t' 2 /3t and U = 9t. 

calculations of the Drude weight for the Hubbard model 
on an isotropic triangular model (c.f., Fig. 4 of Ref. 1251 ). 
The fact that gt < 1 leads to a reduced superfluid stiff- 
ness as is observed in ^LOS. However, if we interpret 
our results in the simplest manner [27| they suggest that 
the most correlated materials have the most strongly sup- 
pressed superfluid stiffness which is the opposite trend to 
that found experimentally 0>EI- 

Ak = J2k' Vk-k'(ck>lC-vi) RV B = 5t^k is the super- 
conducting order parameter. Fig. [21 shows the mean 
of | Ah | and |A k | as functions of U for several values of 
t'/t. A k ^ A k indicates a pseudogap, which is predicted 
to be largest near the Mott transition and have d x 2_ y 2 
symmetry. The angle dependence of the pseudogap could 
be measured by angle resolved photoemission or angle re- 
solved magnetoresistance oscillations. We plot the mean 
of |A k | and |A k | as functions of t'/t for fixed [21 Varying 
t'/t can lead to a suppression of superconductivity and 
can even drive the Mott transition as can be seen from 
the phase diagram of the model (inset to Fig. [2J) . 

Our results suggest that the effects of the anisotropy 
of the triangular lattice are important for the organic 
superconductors. Most importantly we suggest that the 
large value of t'/t and hence of J'/ J stabilises the RVB 
state in ^LOS [24|. The RVB state naturally explains the 
first order transition between the Mott insulator and a 
d x 2_ y 2 superconductor. Further, our results suggest that 
the simple picture 0, 0] in which the only role of hydro- 
static and 'chemical' pressure is to vary U /W is not suf- 



ficient to explain the phase diagram of ^LOS. It appears 
that the value of t'/t, and hence J' /J also plays a cru- 
cial role in determining the behaviour of these materials. 
t'/t controls the degree of nesting of the Fermi surface 
and therefore directly controls the stability of the Mott 
insulator, whereas J' / J determines the pairwise poten- 
tial Vk which controls the stability of the superconduct- 
ing phase. This is why variations in t'/t can even drive 
the Mott transition at a fixed U. Clearly, the physics of 
the anisotropic triangular lattice is qualitatively different 
from that of square lattice. 

Combining our results with those from DMFT stud- 
ies of the Hubbard model allows one to reproduce the 
main features of the phase diagram of ^LOS (Fig. 
This is consistent with the claim that the Hubbard 
model is the minimal model for ^LOS. (We stress that 
this theory is not applicable to the quarter-filled layered 
organic superconductors 01 or the Bechgaard salts 
However, caution is required here. Although calculations 
(e.g., [lfij l based on the Hubbard model can give good 
quantitative agreement with experiments on ^LOS one 
does not know a priori how to map the experimental pa- 
rameter space (pressure, temperature and chemical com- 
position) onto the theoretical parameter space (t'/t, U/t, 
etc.). Therefore an outstanding problem is to discover 
whether quantum chemistry predicts the large variations 
of the Hubbard parameters with pressure required for 
quantitative agreement with experiment. 

We have presented an RVB theory of the Hubbard- 
Heisenberg model on the anisotropic triangular lattice. 
We argued that the RVB state may be a good trial wave- 
function for iLOS because the values of J' /J deduced 
from quantum chemistry and experiment are compara- 
ble to those for which the RVB state appears to be a 
good approximation. Our calculations show a first order 
Mott-Hubbard transition from an insulating state to a 
d x 2_ y 2 superconductor. A similar first order Mott transi- 
tion is seen in experiments on |LOS. The Mott-Hubbard 
transition can be driven by increasing either U /t or t'/t. 
Further, at a fixed U, superconductivity is strongly sup- 
pressed by increasing t'/t. This suggests that the value of 
t'/t may be more important in the layered organic super- 
conductors than has previously been appreciated. The 
superconducting state has a reduced superfluid stiffness 
as is observed in the ^LOS. The RVB theory predicts 
that there is a pseudogap with d x 2_ y i symmetry 

Note added: after completing this work we became 
aware of some similar results obtained by Gan et al. |29| 
and Liu et al. |3(ij |. 
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